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Since the advent of low-lightdevel television {LLLTV) systemsy it has been recognized 
that such devices offer the ability to observe meteors as faint as 10th magnitude iihich 
allows the extension of optical meteor data to masses as small as 10^^ gram. The Space 
Sciences Laboratory at Marshall Space Flight Center has been actively engaged in such 
observations using image orthicons and intensified SEC mdicons, ^ 

The results of these observations are presented along with an interpretahon in terms 
of mass -flu x. This mterpretatioii includes the development of a relationship between 
peak luminosity of a meteor and masSy velocity j and zenith angle that was derived froni 
single body meteor theory and compares favorably with results obtained from the Arti- 
ficial Meteor Program. Also included in the mass-flux inter pr elation is an analysis of 
the observation response of a LLLTV system to fixed and moving point sources. 


T he present model of the meteoroid mass 
distribution is based on an extrapolation 
from ground-based photographic observations of 
the larger meteoroids with masses of the order of 
grams to satellite-borne penetration measure- 
ments of meteoroids with masses in the microgram 
range. The mass range representing the greatest 
damage potential to manned space vehicles is 
from 1 to 100 mg. The fact that the meteoroid 
population in this region must be inferred from 
an extrapolation over 6 orders of magnitude 
between two points that are determined by com- 
pletely different properties of meteoroids through 
interactions that are poorly understood physically 
and cannot be adequately tested experimentally 
has caused some concern among those responsible 
for establishing the meteoroid environment. Also, 
there arc legitimate scientific reasons for ex- 
tending the ground-based optical measurements 
to fainter meteors. Of primary interest is the 


determination of the slope of the mass distribution 
curve or the population index parameter. This will 
greatly improve the confidence in the extrapola- 
tion as well as reduce the range over which the 
extrapolation must be carried out, and will 
provide a badly needed consistency check betw con 
the ground-based and satellite measurements. 

This paper discusses the techniques and results 
of uring LLLTV observations to determine the 
meteoroid mass distribution in the region from 
grams to milligrams. 

DESCRIPTION OF THE SYSTEM 

Recent developments in low-light-level tele- 
vision systems have allowed the observation of 
much fainter meteors than could be photographed. 
This improvement results primarily from the 
much higher quantum efficiency of the photo- 
dctector which results in much smaller integration 
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times. This is particularly important in meteor 
work since it is advantageous to have the integra- 
tion time shorter than the event duration. Our 
system consists of a Commercial Electronics 
camera chain using a Westinghouse WLr32000 
Intensifier-SEC vidicon tube with a 105-mm 
//.75 Rayxar lens. This affords a 13 by 16° field of 
view. The effective integration time of the system 
is very close to the standard frame time of s, 
which is ideal for meteor work. 

The ultimate theoretical sensitivity for the 
system is 14.26. This was estimated by 
requiring a star to produce 1 photoelcctron at the 
photocathode per integration time. The system is 
invariably limited by sky background which, even 
under ideal conditions, is 2 orders of magnitude 
above dark current. Stars as faint as 11 have 
been observed. This is close to the theoretical limit 
of 11.6 for a sky background of 300 Wt =10 
stars/deg2 and a S/N = 5. 

For moving objects, such as meteors, this 
limiting magnitude would apply to meteors 
moving nearly parallel to the optical axis so that 
they remain within a resolving element for one 
integration time. For most meteors, the S/N will 
be decreased because the time they contribute to 
a resolving element is limited by the writing speed 
of their image. It is estimated that if the system 
has an observing limit of 11 magnitude for stars, 
it will see all meteors brighter than w„ = 6.4 and 
50 percent of the meteors brighter than m^=S.lo, 

ANALYSIS PROCEDURES 

Even though the dynamic range of a TV image 
is limited, photometry" of point images can still be 
performed over as much as 6 orders of magnitude 
by making use of the fact that the image spreads 
after it reaches saturation (Beyer et al., 1966). 
Thus, the amount of light associated with the 
image is a monotonic, if not linear, function of 
input. The difficulty lies in the fact that obtaining 
light curves from the TV monitor is a time con- 
suming and laborious task, especially for the faint 
meteors. Until the special video processing systems 
presently being developed for this purpose arc 
available, it will not be possible to obtain light 
curves on a sufficient sample of meteors to 
establish a good distribution. 

An alternative procedure, which is less time 


consuming, was adopted for interim use. This 
consists of treating the video system as a threshold 
detector and simply counting those meteors that 
are above the detection threshold. By varjung the 
threshold through reduction of the lens aperture 
setting, a cumulative distribution in peak meteor 
magnitude is obtained. 

Figure 1 shows the results of this mode of 
operation during two observing periods at Climax, 
Colorado. The camera was oriented toward the 
zenith. Care was taken to program the aperture 
settings to assure a uniform distribution of 
observing time at each aperture setting through- 
out the night. The higher rates in the period from 
July 31 through August 9 may be attributed to 
the presence of the Perseids and the 5-Aquarids 
during this period. The details of these data will 
be published separately. 

Taking the area of the sky wfithin the field of 
view to be 445 km^, the observed rates for the 
non-shower period can be expressed 

log ^ = - 15.352+. 5053m. (1) 

where is the observed rate of meteors (number/ 



Figure 1. — Obsen^cd meteor rates as a function of limiting 
magnitude. Error bars represent la limits based on 
the statistical sample. 
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mVs) when is the limiting magnitude. This 
limiting magnitude is taken to be the magnitude 
of the faintest AO star detectable on the video 
screen. In a sense, this limiting magnitude repre- 
sents the faintest meteor that could be detected, 
although practically all detectable meteors will 
have to be brighter than this because of their 
writing speed. 


RELATIONSHIP BETWEEN PEAK 
LUMINOSITY AND METEOR MASS 


A statistical analysis by Jacchia, Verniani, and 
Brig^ (1965) of the peak brightness" of meteors 
in terms of their mass, velocity, and entry angle 
yielded the result 

7p= (2) 

where /p is the peak intensity in units of zero 
magnitude stars, Wo is the initial meteor mass in 
grams, v is the entry velocity in km/ s, and 6 is the 
entry angle or the angle between the velocity 
vector and zenith. 

A similar result can be derived from classical 
single body meteor theory by assuming that the 
rate of mass loss is equal to the energy input 
divided by the heat of vaporization L, In the free 
molecular regime, 

Apt^ 

^ 2LpJI^ 

where S is the shape factor (jS=L 208 for a 
sphere) , pm is the density of the meteoroid, and 
p is the atmospheric density. Thi.s neglects radia- 
tion losses, which arc small compared to the heat 
input even at the boiling point of Fe (3160° K) ; 
and neglects anj’ change in shape fact or w ith 
time. Neglecting deceleration, which amounts to 
only a few percent velocity change in the time it 
takes a small meteoroid to completely burn up, 
and assuming an exponential atmosphere, equa- 
tion (3) becomes 




2LpJ^ 


pjier^^^ 


( 4 ) 


where — vi cos h is the scale height, and p// 
is the atmospheric density at x = 0. This difT(‘r- 
ential equation may be solved by separation of 
variables. It is convenient to define ^ = 0 at meteor 


burnout, or m = 0. The solution is then 
m = 77io{l — e^*y 

where 


cos 6 



( 5 ) 


and mo, the initial mass, is 



Differentiating the solution, 

771= -307nQ{l-e^^y^^ ( 6 ) 

The peak rh is obtained by equating in from 
equation (6) to 0 to find the time tp when 7h is 
maximized. This yields 

e^'V = H 

Putting this in equation (5) , 

Mp=—%l3mo (7) 

The radiant intensity from a meteor is given by 




T . 9 

- mv- 

2 


( 8 ) 


where r is the luminous efficiency. Using equation 

( 7 ), 


r 4 


COS 6 


( 9 ) 


The luminous efficiency for Fe has been deter- 
mined experimentally from Trailblazer (Ayers 
et al., 1970) and is expressed by 

r= 10*17 *"^Kcm/s) (19) 

Figure 2 compares the light curve obtained with 
this model with an observed light curve from one 
of the Fe Traffilazer meteors. 

For stony meteors. Cook, Jacchia, and 
McCrosky (1963) recommend 

T= 10-1^ -"ly (cm/s) (11) 


which is consistent with estimates of Ayers, 
McCrosky, and Shao based on Trailblazer meas- 
urements. Putting this value in equation (9), and 
choosing A = 5.4 km which corresponds to 80 km, 
results in 


7 p= cos 6 (12) 
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Figure 2. — Comparison of theoretical light curve com- 
pared with the measured light curve of an artificial 
meteor. 


where Ar is the area of a resolving element, F is 
the focal length of the lens, and o> is the angular 
rate. 

The criterion for detectability is 
^ 1/2 

where It is the threshold intensity for stationary 
objects, II is the limiting intensity for moving 
objects, and r is the integration time. The pho- 
tocathode is 40 mm in diameter, or 24 mm from 
the top to bottom raster line. Since there are 525 
lines, a resolving element is taken to be a square 
24/525 mm per side. Taking F=105 mm and 
't = Hq s, 

Il = It 76 . 57 oo; co> 0.013(W rad/s (15) 


vith the units the same as in equation (2). For 
the case of a 1 g meteor at 22 km/s, typical of 
cases from which equation (2) was obtained, 

(1.154 (eq.2) 

/.= ] 

(1.182 (eq. 12) 

The fact that the simple theory yields almost 
identical results as the empirical approach, 
together with the desirability of having some 
theoretical basis for determining the functional 
relationship of luminous intensity with mass, 
velocity, and entry angle, are the bases for 
choosing equation (12) as the functional observing 
relationship for the analysis. 


SYSTEM RESPONSE TO MOVING TARGETS 

As was stated previously, the detection threshold 
of a moving object is increased because the 
photons are spread over a number of resolving 
elements instead of contributing to the signal in a 
single element. A first order attempt to derive the 
system response to moving targets is to simply 
require that the minimum detectable moving 
object deliver the same number of photons during 
the time it resides in a resolving element as a 
minimum detectable stationary object in one 
integration time. The residence time is 



(13) 


Figure 3 compares this model of the response to 
moving point sources with measurements made in 



Figure 3. — Response of the SEC Vidicon to moving point 
images. The solid curves represent the predicted re- 
sponse based on the detection criteria used in the 
analysis. 
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our laboratory. The mcasurementvS consisted of 
projecting point sources onto a mirror mounted 
on a rate table which reflected the image onto a 
projection screen. The point source consisted of a 
precision 50 m pinhole mounted in a 35-mm slide 
projector. Neutral density filters were used to vary 
the image intensity, and the various images were 
calibrated by a photomultiplier photometer. The 
background lighting was provided by a small flood 
lamp with a controllable aperture. The angular 
rate of the mirror was adjusted for each filter until 
the image could no longer be detected on the 
monitor screen. 

The agreement between the simple model and 
the observed results justifies the model. The only 
discrepancy occurs in the transition at the point 
where the image spends r in a resolving element. 
Since the image in the experiment is smaller than 
the resolving element but is not an infinitesimal, 
some rounding off of the theoretical model in this 
transition region is to be expected, which, of 
course, is the observed result. 

Since the angular rate 


T 


and using r = 80 km typical of meteor heights, 
equation (15) becomes 

{ 7r0.957y sin d 

(17) 

It 

whichever is greater. 


RELATION BETWEEN INCIDENT FLUX 
AND OBSERVED FLUX 


only the observed 4>, clearly there is not sufficient 
information to solve the integral equation. 
Several simplifying assumptions are in order. 
First, to a good approximation the velocity and 
mass distributions are independent. Second, it 
will be assumed that the velocities are isotropically 
distributed, i.e., n™. is independent of 6 and <j>. 
Equation (18) becomes 

.1/2 r” f" 

(f> = 2ic f cos 6 wi odd I n,dv J n^dni (19) 

Even assuming 71^ is known, there still is insuffi- 
cient information to define nw However, the fact 
that the observed (f> can be expressed as a power 
law, equation (1), suggests the cumulative mass 
flux Nm (number per unit area time with mass m 
or greater) can be expressed 

N^=tt f Nfndm = Cmr~^ 

'' mT 

where a is the population index. W ith this assump- 
tion equation (19) becomes 

(f, = 2C j cos^sin^d^f n^dimiT^^{lTjOjV) (20) 
do 

The threshold mass is from equation (12) 

But from equation (17) 

e<So 

( 21 ) 

0.957/rl0®*®^i^"^ tan 9 e>do 

where 

(oiw) 


The number of meteors observed per unit area 
time is given by 

<t>= J do cos 6 

hemisphere 

where Umv is the directional mass velocity dis- 
tribution (number per unit area, time, solid angle 
with masses between m and my dm and velocities 
between v and v+dv). The integration is carried 
out from the threshold mr required to produce an 
observed signal, which is a function of the thresh- 
old response of the system 7r, 9j and v. Given 


n 


dv I d?7lTlmv 

mT(7T,SfV) 


(18) 


Equation (20) becomes 

0 = 2C7r-^lO^"*^« 


r y .00 

v^^Hvdv J (cos ^)^+® sin (9 d0 

+ (0.957)"“ f v^’^Uvdv 
do 

/ ir/2 1 

(cos 0) (sin 9) d9 (22) 

-o(») d 
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It is convenient to introduce an average m 
defined as the mass of a just detectable meteor 
having average velocity and 5 = 45°. From equa- 
tion (21) 


7n = 0.957/rl0“ tan 45° 


Equation (22) becomes 

2C7«-“ r 


/ 

^ ft 


dv 


(23) 



(cos 5) ’+“ sin 6 



dv 


•/ (cos 5)*+“(sin 5)‘-“ d5 (24) 

■'*o(») J 

The first integral is the contribution from those 
meteors moving nearly along the line of sight that 
remain in a single resolving clement for one 
integration time. Since do is typically 3°, this con- 
tribution is small and can bo ignored. The integral 
over 8 in the second integral must be evaluated 
numerically unless a is an integer. The lower limit 
5o(i’) is a function of v, however, as may b(! seen in 
figure 4. This dependence is not strong for a«l. 
Therefore, an average value of r = 20km/s will 
be used which yicJds 5o = 3°. 

The velocity distribution was adopted from the 
work of Dohnanyi (1966), exprc.ssed as 


C.vl.61X10V'>3 


11.2<i'<16.6 

(25) 

16.6 <r< 72.2 


The normalization constant 001 153. The 

(v) using this distribution is 19.2 km/s. The 
weighted average (r-')/(t>)^ is shown in figure 5. 
Equation (24) can be written 


<t> = 2 


(w)®“ 


I{6o)N^ 


(26) 


where A^«j. is the cumulative isotropic mass flux of 
meteoroids having mass nir or greater, ?n is given 
by equation (23), and 7(5) is 



Figure 4.— The value of the integral I (Bo) in equation (26) 
as a function of $o and a. 


Since It = 10~”-^ and the observational results 
were found to be as follows, from equation (1) ; 

log ^ = - 15.352 -b0.5053m„ 

log ^=-15.352- 1.263 log 7r (27) 

Differentiating the log of equation (22) 
d(Iog<j>) 

d(log Jr) “ 

From equation (27) a is found to bo 1.263. Using 
this value in figures 4 and 5, 


7(5) =0.933 

Combining equations (23), (26) and (27) 



(cos 5)‘+“(sin 5)'-“ dO 


logAW= -15.352-1.263(logWr-1.466) (28) 

—log 2— log 2.48— log 0.933 
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Figure 5. — Tlic value of the moments of the velocity dis- 
tribution computed from Dohnanyi’s velocity 
distribution. 


or 

log AW= - 14.24-1-263 log mr 

Th(! observed range from ??ip = 7 to mi, = 10 corre- 
sponds to a mass range of to 10'- g. The 

results of equation (28) are compared with the 
existing distribution of meteors in figure G. The 
result from the July 31 to August 9 expedition is 
also shown which contains the Perseids and 
5-Aquarids. No attempt was mad(‘ to alt(T the dis- 
tributions in velocity and angle to account for 
the shower component. 

COMPARISON WITH OTHER DATA 

Figure G shows the relationship of the data 
obtaiiu^d in this work with the current niet(‘oroid 
mass distribution adopted by NASA (Coiir- 



Figure 6.— Comparison of the results of this study with 
other work and with the adopted NASA meteoroid 
environmental design criteria. The observed data 
represent encounter frequency, whereas the design 
criteria are weighted to express penetration frequency. 


Palais, 1969). Also shovTi for comparison is the 
Hawkins and Upton (1958) datum point based 
on photographic meteors, and the points obtained 
from the Pegasus (Clifton and Naumann, 19G6) 
and Explorer XXTII penetration experiments 
(D'Aiutolo, 19G5) . These data have been analyzed 
in terms of encounter frequency, i.c., number 
incident per unit area time on a surface v ithout 
regard to angle of incidence or velocity. The 
NASA design curve refers to the number capable 
of penetrating a surface per unit area time which 
can just be penetrated by a meteoroid with the 
specified mass under conditions of normal impact 
at the average velocity (assumed to be 20 km/s) 
and average density (assumed to be 0.5 g/cm ) . 
Integrating over the velocity and angular dis- 
tributions weighted appropriately for penetration 
mechanics results in the penetration frequency 
being less than the encounter frequency by a 
factor of approximately 2 (Naumann, 1966) . 

The penetration data from Pegasus and Ex- 
plores XXIII were analyzed using recent calibra- 
tion data (Naumann ct ah, 1969). The dis- 
crepancy' bctwi'en these data points and the N ASA 
model is partially due to the difference between 
encounter frequency" and penetration frequency , 
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but also includes the departure from linear size 
scaling in thin metallic targets. Again, the design 
curve was derived from the actual penetration 
data, using the conventional pent^tration formula, 
and will serve adequately as a design criterion so 
long as the same formula is used to convert back 
to penetration results. However, this departure 
sliould be considered in developing a true mass 
distribution. 

CONCLUSIONS 

A technique has been developed using the SEC 
vidicon LLLT\ system as a threshold detector for 


faint meteors to obtain mass flux distribution data 
in the mass range from 1 to 100 milligrams. The 
analysis technique is based on peak intensities 
using the most recent values of luminous efficiency 
obtained from the Trailblazcr measurements. The 
data are quite consistent with present pho- 
tographic data at 1 gram and the satellite data at 
1 microgram, and tend to confirm the adopted 
NASA meteoroid model. 

It is recognized that these observations repre- 
sent only one time during the year and may be 
subject to seasonal variations. Such effects are 
the object of a current investigation. 
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